Search for the Decay jc^uiy 
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We performed a search for the K^ t^^vV decay at the KEK 12-GeV proton synchrotron. No 
candidate events were observed. An upper limit on the branching ratio for the decay was set to be 
6.7 X 10"^ at the 90% confidence level. 

PACS numbers: 13.20.Eb, 11.30.Er, 12.15.Hh 



The decay tt^uu occurs via loop diagrams that 

change the quark flavor from strange to down [1]. It vio- 
lates CP symmetry directly, and the amplitude is propor- 
tional to the imaginary part of the Cabbibo-Kobayashi- 
Maskawa matrix elements in the standard model (SM). 
Since the theoretical uncertainty in the branching ratio is 
small and controlled, it provides a good testing ground of 
the SM and beyond [2]. The branching ratio Br(K^ 
TT^ujy) is predicted to be (2.49 ± 0.39) x IQ-^^ [3]. The 
current experimental limit is 2.1 x 10~^ at the 90% con- 
fidence level by our previous search [4]. 



KEK E391a is the first experiment dedicated to the 
TT^uu decay. Neutral kaons were produced by 
12 GeV protons incident on a 0.8-cm-diameter and 6- 
cm-long platinum target. The proton intensity was typ- 
ically 2 X 10^^ per spill coming every 4 sec. The neutral 
beam [HI, with a solid angle of 12.6 /istr, was defined by 
a series of six sets of collimators and a pair of sweeping 
magnets aligned at a production angle of 4 degrees. A 
7-cm-thick lead block and a 30-cm-thick beryllium block 
were placed between the first and second collimators to 
reduce beam photons and neutrons. The momentum 



peaked around 2 GeV/c at the entrance of the detector, 
11m downstream from the target. 

Figure [T] shows the cross-sectional view of the E391a 
detector, i^^'s entered from the left side, and the detec- 
tor components were cylindrically assembled along the 
beam axis. Most of them were installed inside the vac- 
uum tank to minimize interactions of the particles be- 
fore detection. The electromagnetic calorimeter, labeled 
"Csl" , measured the energy and position of the two pho- 
tons from TT^ . It consisted of 496 blocks of 7 x 7 x 30 cm^ 




FIG. 1: (color online) Schematic cross-sectional view of the 
E391a detector. "Om" in the scale corresponds to the entrance 
of the detector. 
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undoped Csl crystal and 80 specially shaped Csl blocks 
used in the peripheral region, covering a 190 cm-^ cir- 
cular area. To allow beam particles to pass through, 
the calorimeter had a 12 x 12 cm^ hole at the center. 
The main barrel (MB) and front barrel (FB) counters 
consisted of alternating layers of lead and scintillator 
sheets with total thicknesses of 13.5 Xq and 17.5 Xq, 
respectively, and surrounded the decay region. To iden- 
tify charged particles entering the calorimeter, scintil- 
lation counters (CV) hermetically covered the front of 
the calorimeter. It consisted of a plastic scintillator ho- 
doscope that were placed 50 cm upstream of the calorime- 
ter and four 6-mm-thick scintillator plates that were lo- 
cated parallel to the beam axis between the hodoscope 
and the calorimeter. Multiple collar-shaped photon coun- 
ters (CCOO, CC02-07) were placed along the beam axis 
to detect particles escaping in the beam direction. CC02 
was a shashlik type lead-scintillator sandwich counter, 
and was located at the upstream end of the decay re- 
gion. CC03 filled the volume between the beam hole and 
the innermost layers of the Csl blocks in the calorimeter. 
The vacuum region was separated by a thin multi-layer 
film ("membrane") into the beam and detector regions. 
This kept the decay region at 1 x 10~^ Pa despite some 
outgassing from the detector. Further descriptions of the 
E391a detector are given in [3, Ig^. 

The E391a experiment started taking data in Febru- 
ary 2004. In the first period, whose partial analysis was 
reported in 0], the membrane drooped into the neutral 
beam near the calorimeter and caused many neutron- 
induced backgrounds. After fixing this problem, we re- 
sumed the physics run in 2005. In this analysis, we used 
the data in the second period from February to April 
2005. Data were taken with a hardware trigger requir- 
ing two or more shower clusters in the calorimeter with 
> 60 MeV. We also required no activity in the CV and 
in some other photon counters. 

Analysis procedures were basically the same as de- 
scribed in [4]. First, we identified clusters in the calorime- 
ter, each of which should have transverse shower shape 
consistent with a single photon. The clusters were re- 
quired to have more than 150 and 250 MeV for the lower 
and higher energy photons from tt^, respectively. To 
improve the accuracy of energy measurement, clusters 
within the 36 x 36 cm^ square around the beam and 
outside the radius of 88 cm were not used as the pho- 
ton candidates. Second, we selected events with exactly 
two photons in the calorimeter and without any in-time 
hits in the other counters. In order to achieve high effi- 
ciency of particle detection [7.] , energy thresholds for the 
counters were set at around 1 MeV; e.g. 1.0 MeV for FB, 
MB, and CC02, and 0.3 MeV for CV. An additional pho- 
ton counter (BA) was placed in the beam at the down- 
stream end. It consisted of a series of alternating layers 
of lead, quartz, and scintillator plates. Photons to the 
BA were identified by the Cerenkov light in the quartz 



layers and by the energy deposition of more than 20 MeV 
in the scintillator layers. Third, assuming that two pho- 
tons came from a tt^ decay on the beam axis, we calcu- 
lated the decay vertex position along the beam axis (Z) 
and the transverse momentum of tt^ (Pt)- Fourth, we 
imposed kinematic requirements as follows. To remove 
77 decays, we calculated the opening angle be- 
tween two photon directions projected on the calorimeter 
plane, and required it to be < 135 degrees. The shower 
shape was required to be consistent with a photon en- 
tering the calorimeter with the direction from the decay 
vertex to the hit position on the surface. The recon- 
structed TT^ should have the energy less than 2 GeV, and 
should be kinematically consistent with a tt^vv 
decay within the proper momentum range. Finally, 
we defined the region for the candidate events (signal 
box) in the Pt vs Z plot as 0.12 < Pt < 0.24 GeV/c 
and 340 < Z < 500 cm. In this analysis, we masked the 
signal box so that all the selection criteria (cuts) were 
determined without examining the candidate events. 

There were two types of background events. One was 
the events from decays and the other was the events 
due to the neutrons in the halo of the neutral beam ( "halo 
neutrons"). 

The main background source from decays was the 
TT^TT^ mode, whose branching ratio is 8.7 x 10~^. 
There are four photons in the final state, and if two of 
them escape detection, tt^tt^ can fake a signal 

event. The number of background events was estimated 
by Monte Carlo simulation. We generated tt^tt^ 
decays with 11 times larger statistics than the data. Af- 
ter imposing all the cuts, the background level was es- 
timated to be 0.11 events. It turned out that the Csl 
calorimeter and the main barrel were the most responsi- 
ble components for detecting extra photons and rejecting 
those backgrounds. To verify the detection inefficiency of 
photon counters in the simulation, we analyzed the events 
with four photons reconstructed in the calorimeter. As 
shown in Fig. [21 in addition to the tt^tt^ events at 

the mass, there was a tail in the lower mass region due 
to contamination from — > 37r^, whose two out of six 
photons escaped detection. The number of events in the 
tail was reproduced by our simulation [1]. For charged 
decay modes {K^ tt+tt^tt^, nlvil = e,/i)), we 

studied the rejection power of the kinematic cuts to these 
backgrounds. Multiplying the expected inefficiency of 
charged particle counters to their rejection, we estimated 
their contribution to be negligible. 

Halo neutrons induced a substantial portion of back- 
grounds, although the halo was suppressed by 5 orders 
of magnitude from the beam core. The background was 
categorized into three types and they were estimated sep- 
arately. The first type was due to tt^'s produced in the in- 
teraction of halo neutrons with the upstream CC02 collar 
counter ("CC02 BG"). Ideally, their Z position should be 
reconstructed properly at CC02, i.e., outside the signal 
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FIG. 2: Reconstructed invariant mass distribution of the 
events with four photons in the calorimeter. The points show 
the data, and the histograms indicates the contribution of 
7r°7r° and Stt^ decays (and their sum), expected 

from the simulation, normalized with the number of events in 
the Kl 



TT^TT^ peak. 



box. However, they can enter the signal region when the 
energy of either photon was mismeasured due to shower 
leakage or photo- nuclear interactions in the calorimeter. 
To reproduce the tail in the vertex distribution, we used 
data obtained in a dedicated run for the study ( "Al plate 
run"), in which a 0.5-cm-thick aluminum plate was in- 
serted to the beam at 6.5 cm downstream of the rear end 
of CC02 [9]. After imposing analysis cuts and selecting 
events with two photons in the calorimeter whose invari- 
ant mass was consistent with tt^, we obtained the distri- 
bution of reconstructed Z vertex of tt^'s produced at the 
Al plate. It was then convoluted with the Z distribution 
of TT^'s production points within CC02 p^] so as to match 
the peak position with that observed in the physics run, 
as shown in Fig. [3l The distribution was normalized to 
the number of events in Z < 300 cm. We estimated the 
number of CC02 BG events inside the signal box to be 
0.16. 
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FIG. 3: (color online) Reconstructed Z vertex distribution of 
TT^'s produced within CC02. The points show the data in the 
physics run in the upstream region (Z < 340 cm) and the 
histogram indicates the distribution from the Al plate run. 



The second type of neutron-induced background was 
due to neutron interactions with the CV ("CV BG"). 
This background should also be reconstructed properly 
at the Z position of the CV, i.e., outside the signal box. 



However, events can shift upstream when either cluster 
was overlapped by other associated particles and thereby 
mismeasured, or when one of the clusters (or both) was 
in fact not due to a photon from tt^. In order to evaluate 
the background level inside the signal box, we performed 
a bifurcation study with data 0, [l2[. In the simula- 
tion studied beforehand, the cuts against extra particles 
and the shower shape cut turned out to be efficient in 
the background reduction; these cuts were chosen as two 
uncorrelated cut sets in the bifurcation study. The re- 
jection power of one cut set was evaluated with inverting 
another cut set, and vice versa. Multiplying the obtained 
rejection factors, the number of CV BG events inside the 
signal box was estimated to be 0.08. 

The third type of neutron-induced background was due 
to 77's produced by the halo through interactions with the 
CV ("CV-7^ BG"). Since the Z vertex position was cal- 
culated by assuming the tt^ mass, r^'s were reconstructed 
about four times farther away from the calorimeter, and 
they can fall into the signal box. To simulate the rj pro- 
duction, we used a GEANT4-based simulation with the 
Binary Cascade hadron interaction model [13]. Figured 
demonstrates the simulation, which reproduced the in- 
variant mass distribution (from tt^ mass to rj mass) of 
the events with two photons in the calorimeter from the 
Al plate run, normalized by the number of protons on 
the target. We then simulated rj production at the CV 
and estimated the number of CV-r^ BG events inside the 
signal box to be 0.06. 
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FIG. 4: (color online) Reconstructed invariant mass distribu- 
tion of the two photon events in the Al plate run, explained 
in the text. Points with error bars show the data. Histograms 
indicate the contributions from 7r° and 77 produced in the Al 
plate, 77 decays, and their sum, respectively, from the 

simulation. Events in the low mass region were considered to 
be due to neutron interactions accompanying neither tt^'s nor 
?7's, which were not recorded in the simulation. 



Table [J summarizes the estimated numbers of back- 
ground events inside the signal box. We also examined 
the numbers of events observed in several regions around 
the signal box, and they were statistically consistent with 
the estimates. 

After determining all the selection criteria and esti- 
mating background levels, we examined the events in the 
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TABLE I: Estimated numbers of background events (BG) in- 
side the signal box. 



Background source 


Estimated number of BG 


K]^ — > 71 71 


0.11 ±0.09 


CC02 


0.16 ±0.05 


cv 


0.08 ±0.04 


CV-ry 


0.06 ±0.02 



total 0.41 ±0.11 
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signal box and found no candidates, as shown in Fig. O 
The number of collected decays was estimated us- 
ing the TT^TT^ decay, based on 1495 reconstructed 
events, and was cross-checked by measuring Stt^ 
and 77 decays [14]. The 5% discrepancy observed 
between these modes was accounted for as an additional 
systematic uncertainty. The single event sensitivity for 
the TT^vu branching ratio is given by 



Acceptance • A^(i^£decays) 

where the acceptance includes the geometrical accep- 
tance, the analysis efficiency, and the acceptance loss due 
to accidental hits. Using the total acceptance of 0.67% 
and the number of decays of 5.1 x 10^, the single 
event sensitivity was (2.9 ± 0.3) x 10~^, where the error 
includes both statistical and systematic uncertainties. 

Since we observed no events inside the signal box, we 
set an upper limit for the ir^vv branching ratio, 

Br{Kl TT^uu) < 6.7 x 10"^ (90%aL.) , 

based on the Poisson statistics. In deriving the limit, 
the uncertainty of the single event sensitivity was not 
taken into consideration. The result improves the previ- 
ous limit [4] by a factor 3, and the background level by 
an order of magnitude. 
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FIG. 5: (color online) Scatter plot of Pt vs reconstructed Z 
position after imposing all the cuts. The points show the data 
and the contour represents the simulated distribution of the 
signal. The rectangle indicates the signal region. 
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